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Abstract 

The potential of coming experiments on Large Hadron Collider (LHC) to observe the 
rescattering and energy loss of heavy quarks in the dense matter created in heavy ion 
collisions is discussed. We analyze the sensitivity of high-mass /U + /x~ pairs from BB 
semileptonic decays and secondary J/^'s from single B decays to the medium-induced 
bottom quark energy loss. 







1 Introduction 



The significant progress in lattice QCD calculations, in particular, including dynamical quarks, 
strongly suggests that the deconfmement of hadronic matter and chiral symmetry restoration 
must happen at temperatures above T c ~ 200 MeV U. Experimental studies of the properties 
of strongly interacting matter at high enough energy densities for a relatively long-lived quark- 
gluon plasma (QGP) to be formed is one of the goals of modern high energy physics (see, for 
example, reviews in Ref. @, §, |], ||)- It is expected that the quark-hadron phase transition, 
which likely occurred during the first few microseconds of the evolution of the universe, can be 
attained in heavy-ion accelerators. 

In recent years, a great deal of attention has been devoted to "hard" probes of the QGP 
including heavy quarkonia, hard hadrons, jets, and high mass dimuons which are not part 
of the thermalized system, thus carrying information about the early stages of the evolution. 
In particular, the predicted charmonium suppression by screening of bound cc pairs ("colour 
dipole" ) in a plasma or dynamical dissociation by semi-hard deconfined gluons is one of 
the most promising signals of QGP formation. A similar phenomenon has been observed in the 
most central Pb+Pb collisions at CERN-SPS ||: the anomalously small J /if) to Drell-Yan ratio, 
inconsistent with pre-resonance absorption in cold nuclear matter. Although the interpretation 
of this phenomenon as a result of QGP formation is plausible, alternative explanations such as 
rescattering with comoving hadrons can not be fully dismissed [g]. For the heavier bb bound 
states (the T family), a similar suppression effect in QGP is expected at the higher temperatures 
which can be attained in heavy ion collisions at the LHC. 

Along with quarkonium suppression one process of interest is the passage of coloured jets 
through dense matter. Dijets are created at the very beginning of the collision process, < 0.01 
fm/c, by the initial hard parton-parton scatterings. These hard partons pass through the dense 
matter formed by minijet production at longer time scales, ~ 0.1 fm/c, and interact strongly 
with the constituents of the medium. The inclusive cross section for hard jet production is 
negligible at the SPS but increases rapidly with the collision energy. Thus these jets will play 
an important role at RHIC (with Au+Au collisions at = 200^4 GeV) and LHC (with Pb+Pb 
collisions at \/s = 5.5 A TeV). The challenge is to determine the behaviour of coloured jets in 
dense matter [|H]] due to coherent medium-induced gluon radiation JTI], |TJ, [T3|, [TJ], [15], [L6| and 
collisional energy loss due to elastic rescatterings |T7|, [18], |19j. Since the jet rescattering intensity 
strongly increases with temperature, formation of a "hot" QGP at initial temperatures up to 



To ~ 1 GeV at LHC |20| should result in much larger parton energy loss compared to "cold" 
nuclear matter or a hadronic gas. 

The following signals of medium-induced energy loss have been identified as being observable 
in ultrarelativistic heavy ion collisions: 
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1) The suppression of high-p^ jet pairs pi] , p2| with a corresponding enhancement of the 
monojet-to-dijet ratio |23|, p4|. Jets are produced in the initial scattering processes such as 



99 ~> 99 , 19 ^Q9 , QQ -> 99 , 00 -> 99 
where the — > process is dominant. 

2) The pr-imbalance between a produced jet with a gauge boson in 7+jet and Z+jet p6 | 
production, e.g. 

Q9 -> 97 , 90 -> 9^- 
The Z-boson is identified through its decays to muon pairs, Z — > 

3) The modification of the high mass dimuon spectra from semileptonic B and D meson decays 



due to bottom and charm quark energy loss [27, 28, 2£ 



gg -> 66 (cc) -> 55 (DD) -> p + /i"X 



The above phenomena can be experimentally studied in heavy ion collisions |30| with Com- 
pact Muon Solenoid (CMS), which is the general purpose detector designed to run at the 
LHC In our previous work |2^] we discussed the impact parameter distribution in the 

hard jet production processes (jet + jet, 7 + jet and Z + jet channels) at LHC, which can 
shed a light on the dependence of parton energy losses on distance traversed (the coherence 
pattern of the medium-induced radiation can result in non-trivial behaviour of such depen- 



dence, see [13]). The main goal of this paper is to analyze the sensitivity of both high-mass 
pairs from BB semileptonic decays and secondary J/^'s from single B decays to 6-quark 
energy loss in CMS conditions. Note that the ALICE experiment |J2J will also study dilepton 
production in heavy ion collisions although with a different rapidity acceptance: \r]^\ < 2.4 
for CMS while |?7 e | < 0.9 and 2.5 < \r]^\ < 4 for ALICE. We believe that dimuon production 
in combination with high-p^ jet production by gluon and light quark fragmentation can give 
important information about the medium-induced effects for both light and heavy partons in 
heavy ion collisions at the LHC. 

The high-mass dimuon spectra in Pb+Pb collisions at LHC has already been estimated 



in the CMS acceptance |29] assuming a constant energy loss per unit length of dE/dx = 1 
GeV/fm. In our paper we consider the dynamical evolution of heavy quark energy loss and 
rescattering as a function of energy density in an expanding "hot" gluon-dominated plasma. 
We also investigate, for the first time, these effects on secondary J/ip production from single 
-B-meson decays. 

The outline of the paper is as follows. In Sect. 2 we discussed the problem of medium-induced 
gluon radiation of massive quark for incoherent and coherent cases. A model for rescattering of 
a heavy quark in dense QCD-medium is described in Sect. 3. Sect. 3 and Sect. 4 are dedicated to 
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the analysis of calculation results for high- mass dimuon and secondary J/ty spectra respectively. 
The summary can be found in Sect. 5. 



2 Medium- induced gluon radiation from massive quarks 



We first discuss the effects of the medium on gluon radiation from heavy quarks. We emphasize 
that the recent theoretical developments on medium-induced gluon radiation |13| , 



141 are valid 



in the coherent domain only in the ultrarelativistic limit of the quark momenta. Although some 
attempts have been made to calculate medium-induced heavy quark energy loss for quarks of 



mass M q ||33[] , a full description of the coherent radiation from a massive object still does not 
exist. There are two extreme limits for energy loss by gluon radiation. In the low pt limit, ^ 
M q , medium- induced radiation should be suppressed by the mass, while the ultrarelativistic 
limit, pr 3> M q , corresponds to the radiation spectrum of massless quarks. 



The approach developed in Ref. |33| is based on a factorization of the matrix elements into 



elastic scattering and gluon emission where the multiplicity distribution of the radiated 
gluons can be written as 



3a,. 



\rj\ 



dr]d 2 q T ir 2 q \[^ T _ l T f 

where q = (qo,q~T,q3) and I = (lo,lT,h) are the 4-momenta of the emitted and the exchanged 
gluon respectively and r\ = (1/2) In [(g + Q3)/(<lo — Q3)] is the rapidity of emitted gluon relative 
to the initial quark momentum. It is valid in a limited gr-region for small rapidities rj ~ 
and lT{qo/E) <C qx- The coherent LPM effect is taken into account simply by including a 
formation time restriction via a step function [121 • However, for light quarks, it leads to a 



different dependence of the radiative energy loss on the initial quark energy E than found in 



other, more recent works |[T3|, |14 



In our case, the main contribution to high-mass dimuon and secondary J/ip production is 
due to 6-quarks with "intermediate" values of pr > 5 GeV/c, expected to be rather close to the 
incoherent regime. In order to estimate the sensitivity of the dimuon spectra to medium-induced 
effects, we consider two extreme cases: (i) the "minimum" effect with collisional energy loss only 
and (ii) the "maximum" effect with collisional and radiative energy loss in the incoherent limit of 
independent emissions without taking into account the LPM coherent suppression of radiation 
(i.e. dE/dx oc E and is independent of path length, L). In the latter scenario we use the 
Bethe-Heitler cross section obtained in relativistic kinematics and derive the medium-induced 



radiative energy loss per unit length fllj] as the integral over the gluon radiation spectrum 



dE 

dx 



1-M q /E 



Ep 



dy 



4a s C 3 (y)(4-4y + 2y 2 



9iry M 2 y 2 + m 2 (l~y) 



(2) 
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= 9,a 2 s C ab r + _ y2] ^ 2tfpEy(l-y))' 

where m 9 ~ 3T is the effective mass of the emitted gluon at temperature T, y = qo/E is 
the fraction of the initial quark energy carried by the emitted gluon, p oc T 3 is the density of 
the medium and the Debye screening mass squared, fi 2 D = Aiia s T 2 {l + Nf/6), regularizes the 
integrated parton rescattering cross section. 



3 Heavy quark rescattering in a dense QCD-medium 

We have developed a Monte-Carlo simulation of the mean free path of heavy quarks (M& = 
5 GeV and M c = 1.5 GeV) in an expanding QGP formed in the nuclear overlap zone in 
Pb+Pb collisions. The details of the geometrical model of hard quark production and the 



quark passage through dense matter can be found in our recent work p2| . In general, the 
intensity of rescattering and energy loss are sensitive to the initial conditions (energy density 
6q and formation time To) and space-time evolution of the medium, treated as a longitudinally 
expanding quark-gluon fluid. The partons are produced on a hyper-surface of equal proper 



times r = \Jt 2 — z 2 |35] . 

If the mean free path of a hard parton is larger than the screening radius in the QCD- 
medium, A ^> /i^ 1 , the successive scatterings can be treated as independent fI2f . The transverse 
distance between successive scatterings, Axj = (r i+1 — t^vt = — r^pr/E, is generated 
according to the probability density 

dP 1 . f l ds 

ex P(" / uZ I J ( 3 ) 



d(Axi) A(t; + i) J A(r< + s 

where the mean free path is A = l/(erp). The density of the medium p(r) and the quark 
rescattering cross section er(r) are functions of proper time. Then the basic kinetic integral 
equation for the energy loss AE as a function of initial energy E and path length L has the 
form 

AE(L,E)=hx d -^X(x) dE ^ E \ (4) 
J dx dx 
o 

dP(x) 1 

— — = — - exp (-x/X(x)) , 
dx a(x) 

where the current transverse coordinate of a quark, x(t), is determined from dx/dr = vt with 
x — r at vt — 1. 

The dominant contribution to the differential cross section da/dt for scattering of a quark 
with energy E and momentum p = J E 2 — M 2 off the "thermal" partons with energy (or 
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effective mass) mo(r) ~ 3T(r) <C E at temperature T can be written in the target frame 
as 0,|3| 

da ab 2ira 2 (t) E 2 

where C a b = 9/4, 1, and 4/9 for gg, gq and qq scatterings respectively. The strong coupling 
constant is 

12tt 

as[t} ~ (33-2iV»ln(t/A| CD ) (6) 
for Nf active quark flavours and QCD scale parameter Aqcd on the order of the critical temper- 
ature, Aqcd — T c . The integrated parton scattering cross section is regularized by the Debye 
screening mass squared p 2 D : 

^max , 

"-w = I dt ^r < 7 > 

where t max = [s — (M q + m ) 2 }[s — (M q — m ) 2 }/s and s = 2m E + m 2 , + M 2 . 

In the i-th rescattering off a comoving medium constituent (i.e. with the same longitudinal 
rapidity Y) with squared momentum transfer t% and effective mass mo*, the quark loses total 
energy Aej and transverse energy Ae^i as well as getting a transverse momentum kick Ak t j 
relative to the initial momentum pt- It is straightforward to evaluate Ae^, Ae-pi and Ak^: 

Ae Tl ~ (8) 
2m 0i 

Aej ~ Ae^j cosh Y" — Afc t j cos sinh Y, (9) 
Ak,, 



^ 2m J V T Pt 2m 0i 2p T ) Mq ^°' ) 

the angle between the direction of vector k t i and axis z being distributed uniformly. The 
medium-induced radiative energy loss is calculated with Eq. (|) without modification of the 
longitudinal rapidity. 



In our calculations, we use the Bjorken scaling solution for the space-time evolution of 
the energy density, temperature and density of the plasma: 

s(r)r^ = , r 4/3 , (11) 
T(t)tV 3 = T r 1/3 , (12) 
p(t)t = p r . (13) 

To be specific, we use the initial conditions for a gluon-dominated plasma expected for central 
Pb+Pb collisions at LHC @: r ~ 0.1 fm/c, T ~ 1 GeV, N f « 0, p g « 1.95T 3 . It is 



interesting that the initial energy density, £o> i n the dense zone depends on b very slightly, 
feo/^o ~ 10%, up to 6 ~ Ra and decreases rapidly for b > Ra p2f . On the other hand, the 
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proper time of a jet to escape the dense zone averaged over all possible jet production vertices, 
(tl), is found to decrease almost linearly with increasing impact parameter. This means that 
for impact parameters b < Ra, where ~ 60% of the heavy quark pairs are produced pD[ , the 
difference in rescattering intensity and the corresponding energy loss is determined mainly by 
the different path lengths rather than the initial energy density. 

The simulation of quark rescattering is halted if one of the following three conditions is 
fulfilled: 

1) A quark escapes from the dense zone, i.e. its path length becomes greater than the effective 
transverse spread of the matter from the production vertex to the escape point. The details of 
the geometrical calculations of these quantities at a given impact parameter can be found in 
Ref . |2 . 

2) The plasma cools down to T c = 200 MeV. We thus neglect possible additional small contri- 
butions to the total energy loss due to re-interactions in the hadron gas. 

3) A quark loses so much energy that its transverse momentum drops below the average 
transverse momentum of the "thermal" constituents of the medium. In this case, such a quark 
is considered to be "thermalized" and its momentum in the rest frame of the fluid is generated 
from the random "thermal" distribution, dN/d 3 p oc exp (-E/T), boosted to the center-of-mass 
of the nucleus-nucleus collision |28|, . 



4 High-mass dimuon production at LHC 

Let us first consider dimuon production in the high invariant mass region, 20 < M M + M - < 50 
GeV/c 2 , where 

= v /(£ M+ + J EV) 2 -(p^+p / ,-) 2 . 

One of the main dimuon sources in this "resonance-free" mass region is the semileptonic decays 
of open bottom and charm mesons. Heavy quark pairs are produced at the very beginning of 
the nuclear collisions by hard gluon-gluon scatterings and propagate through the dense medium. 
They finally form B and D mesons by "capturing" u, d or s quarks during the hadronization 
stage. These mesons will decay with the average meson lifetimes ct b ± = 496 fim, cr^o = 464 
/im, ctd± = 315 /im and cr^o = 124 /xm. We note that ~ 20% of B mesons and ~ 12% of 
D mesons decay to muons. About half of the muons from B decays are produced through an 
intermediate D |36j and contribute to the softer part of the pr-spectrum. There is also dimuon 
production from single B decays: B — > Dfi + X — > [i + fi~Y. The branching ratio for this channel 
is comparable to the yield from bb pair decay. However the muon pairs from single B's are 
concentrated in the low-mass region, M A1 + At - < Mb = 5.3 GeV/c 2 , below our interest here. 
Note that at LHC energies, there can be a significant contribution to heavy flavour produc- 
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tion from gluon splittings, g — > QQ, in initial- or final-state shower evolution [p^j . Although the 
probability for a high-p-r event to contain at least one bb or cc pair is fairly large, most of these 
quarks are carrying a small fraction of the total transverse momentum of the jet and dimuons 
produced in such a way are expected to be concentrated mainly in the low invariant mass re- 
gion. Nevertheless the additional contribution of " showering" heavy quarks can be important 
even in high dimuon mass domain, and this phenomenon is studied in our next work |38|] . 

The main correlated background is Drell-Yan production, qq — > p + pT . The uncorrelated 
part of the dimuon background, random decays of pions and kaons and muon pairs of mixed 
origin, is comparable with the signal from 66-decays |30| but these random decays also appear 



in the like-sign dimuon mass spectra. Thus such background can be estimated from the fi + n + 
and event samples as 



dM V dM dM ' y ' 

and subtracted from the total p + fi~ distribution. 

The cross sections, crjyjy, for heavy quark production in NN collisions at ^/s — 5.5 TeV, the 
initial QQ momentum spectra, and the B and D meson fragmentation have all been obtained 
using PYTHIA5.7 with the default CTEQ2L parton distribution functions and including 
initial and final state radiation in vacuum which effectively simulates higher-order contributions 
to heavy quark production. The corresponding Pb+Pb cross section is obtained by multiplying 
a< NN by t ne number of binary nucleon-nucleon sub-collisions. The initial distribution of QQ 



pairs over impact parameter b can be written as 40 



^2^-0 ,2 AA 

r!0 (b, v G)=T A A(b)a%%(v r s)^^(b,V~s) (15) 



where the differential inelastic AA cross section is 

d 2 a AA 



m 



d?b 



1 



1 - —T AA (b)^ N (V~s) 



(16) 



and the total inelastic non-diffractive nucleon-nucleon cross section is cr^y ~ 60 mb at y/s — 5.5 
TeV. The standard Wood-Saxon nuclear overlap function is Taa(P) = J d 2 sT^^T^db — s|) 
where T A (s) = A J dzpA^s, z) is the nuclear thickness function with nucleon density distributions 
Pa(s,z) [0. 

We also take into account the modification of the nucleon structure functions due to the 
initial state nuclear interactions such as gluon depletion, nuclear shadowing, using the phe- 
nomenological parameterization of Eskola et al. [Q. Figure 1 shows the initial invariant 
mass spectra from correlated bb, cc and Drell-Yan production in the CMS acceptance, > 5 
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GeV/c and \r]^\ < 2.4. The upper and lower histograms for each contribution show the re- 
sults without and with nuclear shadowing respectively. In this kinematical region the influence 
of nuclear shadowing on heavy quarks is relatively small, giving an pa 15% reduction in the 
open bottom and charm decays. Drell-Yan production is somewhat more affected, with an 
~ 25% reduction, because quark shadowing is stronger than gluon shadowing at low x in the 
EKS model [42|. The total impact-parameter integrated rates are normalized to the expected 
number of Pb+Pb events during a two week LHC run, R = 1.2 x 10 6 s, assuming luminosity 
L = 10 27 cm-V 1 to that 

N{^^)=Ra^~L . 

We see that the bb — > rates are greater than those from the other sources by a factor 

of at least 5. Thus we only consider dimuons from bb decays in the remainder of the discussion. 
Moreover, the medium-induced charm quark energy loss can be significantly larger than the 
6-quark loss due to the mass difference, resulting in an additional suppression of the cc — * / u + 
yield. 

Drell-Yan pairs are unaffected by medium-induced final state interactions. These dimuons 
are directly from the primary nuclear interaction vertex while the dimuons from B and D meson 
decays appear at secondary vertices some distance from the primary vertex. The path length 
between the primary vertex and secondary vertices are determined by the lifetime and 7-factor 
of the mesons. This fact allows one cut to suppress the Drell-Yan rate by up to two orders 
of magnitude using the dimuon reconstruction algorithm based on the tracker information on 



secondary vertex position j^Tj . 

Figure 2 presents the invariant mass spectra from bb decays without and with the 

medium-induced 6-quark energy loss described in sections 2 and 3, including nuclear shadowing. 
The initial dimuon rate in the mass range 20 < M„+„- < 50 GeV/c 2 is 2.8 x 10 4 events per 
two week run without energy loss or shadowing. The rate can be reduced by factor of 1.6 — 4 
due to rescattering and 6-quark energy loss in the QGP. Note that the rate integrated over all 
phase space is always conserved: the suppression in the rate appears only when kinematic cuts 
are specified. The dimuon suppression due to collisional loss is more pronounced at relatively 
lower invariant masses since the collisional energy loss is almost independent of the initial quark 
energy. The relative contribution of medium-induced gluon radiation to total energy loss grows 
with increasing M^^- because the radiative energy loss grows with increasing initial energy, 
oc E in the incoherent limit used in our calculations, and thus contributes to the whole high- 
mass dimuon range, especially in the large invariant mass domain. The loss increases for high 
momentum quarks which contribute to the larger dimuon masses. 

Let us mark that the obtained value of dimuon rate suppression factor for the scenario 
with collisional and radiative energy loss (~ 4 at M^+^- ~ 20 GeV) is close to the result 
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of paper [29|, where a constant energy loss per unit length of dE/dx = 1 GeV/fm has been 
assumed. However, there exist the qualitative difference between invariant mass dependence of 
the effect in two models. In the case of constant (independent of energy) loss of heavy quark the 
suppression factor in decreases with increasing dimuon mass (this is similar to our scenario 
with collisional loss only), while the M-dependence in our model is specified by the dominant 
mechanism of energy loss. 

To conclude this section, we note that there are theoretical uncertainties in the bottom 
and charm production cross sections in nucleon-nucleon collisions at LHC energies which affect 
the quark production rates without any loss effects. The absolute dimuon rates depend on 
the parton distribution functions, the heavy quark mass, the 5-meson fragmentation scheme, 
next-to-leading order corrections, etc. Another possible contribution to the b quark cross sec- 



tion could be light gluinos decaying into b squarks P3[ , as proposed to explain the apparent 
deficit in the B cross section at the Tevatron. It is therefore desirable that high mass dimuon 
measurements in pp or dd collisions are made at the same or similar energy per nucleon as in 
the heavy ion runs. 



5 Dimuons from B —> J/vp at LHC 

We now consider another process which can also carry information about medium-induced 
rescattering and 6-quark energy loss: secondary J /if) production. The branching ratio B — > 
J/ipX is 1.15%. The J/ip J s subsequently decay to dimuons with a 5.9% branching ratio |36| so 
that e.g. 

gg ^bb -> BB X -> J/ipY -> ^^Y. 

Figure 3 shows the transverse momentum and pseudo-rapidity distributions of secondary 
J / ip decays for the same nuclear shadowing and energy loss scenarios described in the previous 
section. The total rates are again normalized to the expected number of events in a two week 
Pb+Pb run. We expect 1.3 x 10 4 dimuons from secondary J/ip's at M^+^- = Mj/^ = 3.1 
GeV/c 2 . Including nuclear shadowing reduces this J/ip yield by ~ 25% while the final state 
rescattering and energy loss by 6-quarks can further reduce the J/ip rates by a factor of 1.3 — 2.2 
in the CMS acceptance. 

We see that the influence of nuclear effects on secondary J/ip production and high-mass 
dimuon rates are quite different. The 25% decrease of the secondary J/ip rate by nuclear 
shadowing is comparable to the ~ 30 — 50% effect of medium-induced final state interactions. 
On the other hand, the high-mass dimuon rates are reduced by up to a factor of 4 due to 
energy loss, much larger than the 15% nuclear shadowing correction. The increased sensitivity 
to nuclear shadowing is due to the different x and Q 2 regions probed. The different influence 
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of energy loss on secondary J/ip and high-mass dimuons is because secondary J/ip's come from 
the decay of a single b quark instead of a bb pair and there is a non-negligible probability that 
the energy lost by one quark is small. Thus a comparison between high-mass dimuon and 
secondary J/ip production could clarify the nature of energy loss. 

We note that about 5000 primary <//?/>(—> fi + fi~) in a two week Pb+Pb run are expected 
to be initially produced by gluon-gluon fusion at CMS acceptance. The final primary J/ip 
rate is rather uncertain: on one hand, it should be suppressed due to colour screening || 
and/or dynamical dissociation in a QGP 0; on the other hand, "thermal" models predict some 
additional yield of J/^'s from a QGP, see, e.g. Refs. and references therein. Other models 
suggest that J/ip 's can be regenerated in the hadron phase by DD interactions [fHf) . Note that 
we call all these J/ip's "primary" in the sense that they are from the primary nuclear interaction 
vertex. Although they could be produced late in the collision, either after thermalization or 
in the hadronization stages, the time scale of their formation is much less than the formation 
time of secondary J/^'s from B meson decays. Thus it is necessary to distinguish secondary 
from primary J/^'s. As in the case of separating dimuons from Drell-Yan production and bb 
decays, primary J/i/j's produced at the nuclear interaction point can be rejected using tracker 
information on the secondary vertex position. In addition, primary J/^'s have a softer dimuon 
Pt spectrum than secondary J/iji's and the primary J/ip contribution disappears rapidly with 
increasing px- 

6 Conclusions 

To summarize, we have analyzed the sensitivity of spectra of high-mass n + fi~ pairs from bb 
decays and dimuons from secondary J/ip's to medium-induced bottom quark energy loss in 
Pb+Pb collisions at LHC For certainty the CMS experiment kinematical acceptance has been 
considered. Since a complete description of coherent radiation by massive quarks is still lack- 
ing, the dimuon spectra were calculated for two extreme cases: with collisional loss only, the 
"minimum" effect, and with collisional and radiative loss estimated in the incoherent limit, the 
"maximum" effect. 

We have found that medium-induced parton rescattering and energy loss can reduce the 
dimuon rate in the invariant mass range 20 < M M + A1 - < 50 GeV/c 2 by a factor from 1.5 to 
4 while nuclear shadowing is only of order of 15% in our kinematical region. The relative 
contribution of radiative energy loss to the total dimuon suppression grows with increasing 
invariant mass and pt due to the stronger energy dependence of radiative loss relative to 
collisional loss. 

Since secondary J/ip production reflects the medium- induced energy loss of only one b- 
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quark, the corresponding suppression by a factor of 1.3 — 2 is less than for bb decays. On the 
other hand, the influence of nuclear shadowing on the relatively low invariant mass region, 
M^+ft- ~ Mj/^, seems to be non-negligible. We suggest that comparing high-mass dimuons 
with secondary J/ip production would help clarify the nature of this phenomenon. 

The experimental recognition of high-mass dimuons from bb decays relative to Drell-Yan 
pairs, as well as the recognition of secondary compared to primary J/^'s, could be performed 
using tracker information on the secondary vertex position. 

We conclude that the dimuon spectra will be sensitive to final state rescattering and energy 
loss of bottom quarks in dense matter. However, there are still theoretical uncertainties in 
the initial production of heavy flavours in nucleon-nucleon collisions at LHC energies. Thus 
measurements in pp or dd collisions at the same or similar energies per nucleon as in the heavy 
ion runs are required. 
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Figure 1: The initial invariant mass distribution of pairs from correlated sources with 

Pt > 5 GeV/c and |r/ M | < 2.4. The results are for bb (solid and dot-dashed histograms - without 
and with nuclear shadowing respectively), cc (dashed and dotted histograms), and Drell-Yan 
production (closed and open circles.) 
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Figure 2: Invariant mass distribution of pairs from bb decays with pi^ > 5 GeV/c and 

| r/^l < 2.4 for various scenarios: without energy loss (solid histogram), with collisional loss only 
(dashed histogram), with collisional and radiative loss in incoherent limit (dotted histogram). 
Nuclear shadowing has been included. 
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Figure 3: Transverse momentum and pseudo-rapidity distributions of secondary J ftp decays 
with Pt > 5 GeV/c and |r^| < 2.4 for various scenarios: without nuclear shadowing and 
energy loss (dash-dotted histograms), with shadowing and without loss (solid histogram), with 
shadowing and collisional loss (dashed histogram) and with shadowing, collisional and radiative 
loss in the incoherent limit (dotted histogram). 
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